Consensus of self-driven agents with avoidance of collisions 



o 
o 

m 
(N 



Liqian Peng^, Yang Zhao^, Baomei Tian^, Jue Zhang^, Bing-Hong Wang^'^, Hai-Tao Zhang'^''*, and Tao Zhou^'^^^ 

^Department of Modern Physics and Nonlinear Science Center, 
University of Science and Technology of China, Hefei 230026, P. R. China 
^Research Center for Complex System Science, University of 
Shanghai for Science and Technology, Shanghai 200093, P. R. China 
^ Department of Control Science and Technology, 
Huazhong University of Science and Technology, Wuhan 430077, P.R. China 
* Department of Engineering, University of Cambridge, Cambridge CB2 IPZ, U.K. 
^Department of Physics, University of Fribourg, Chemin du Muse 3, CH-1700 Fribourg, Switzerland 

(Dated: January 23, 2009) 

In recent years, many efforts have been addressed on collision avoidance of collectively moving 
agents. In this paper, we propose a modified version of the Vicsek model with adaptive speed, which 
can guarantee the absence of collisions. However, this strategy leads to an aggregated state with 
slowly moving agents. We therefore further introduce a certain repulsion, which results in both 
faster consensus and longer safe distance among agents, and thus provides a powerful mechanism 
for collective motions in biological and technological multi-agent systems. 
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I. INTRODUCTION 



One of the most marvelous and ubiquitous phenomena 
in nature is collective motion, a kind of motion that can 
be observed at almost every scale: from bird flocks and 
fish schools at a macroscopic level to bacteria, individual 
cells and even molecular motors at a microscopic level 
P, 0, S, IS S S H ■ Ahhough in most cases agents do 
not share global information and often travel in the ab- 
sence of leaders or external forces, collective motion may 
still occur. Analogous behaviors are reported in engineer- 
ing systems also, such as groups of autonomous mobile 
robots and air vehicles [UlllEIilQEIill (see also 
a newly reported swarm model that may connect granu- 
lar materials and agent-based models [l3|)- In order to 
uncover the underlying mechanism leading to the consen- 
sus of collective population, Vicsek et al. proposed 
a model with self-driven agents to mimic the biological 
swarm, which displays a novel type of kinetic phase tran- 
sition. From then on, the nature of the nonequilibrium 
phase transition of collective motion attracted more and 
more attentions [H, [13, [SI, [H [H [Hi . Due to simplic- 
ity and efficiency, many modified versions of the Vicsek 
model were proposed. For example, some new methods 
with effective leadership were introduced [H, [2^, [1^ , and 
new moving protocols with adaptive speed to accelerate 
the consensus were designed [13, [11]; meanwhile some 
scholars have studied the consensus of collective motions 
via low-cost communication 29] and predictive mecha- 
nism [S^, [si], J all of which can greatly enhance the 
global convergence. 
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FIG. 1: Illustration of the current model with adaptive speed, 
where r denotes the horizon radius of an agent, di denotes 
the distance of the ith agent and its nearest neighbor, Vmax 
denotes the possibly maximal velocity, a denotes the size of 
an agent. Accordingly, 2a corresponds to the least distance 
of two agents. 



Recently, much attention has focused on how to keep 
distances among agents. A common way is to introduce 
attraction and/or repulsion [H, [H, [l3, [M [M [IS [H, [13, 
[H, [13, \^ . However, any kinds of repulsions alone can 
not sufficiently avoid collision at all times, because it is 
entirely possible that in a high-density area, two agents 
are compelled to collide for the purpose of avoidance of 
collision with a third agent. In this article, we propose 
a swarm model with adaptive speed to completely elimi- 
nate collisions. In a plane, each agent adjusts its direction 
as the average direction of its neighbors while resets its 
speed according to the minimal distance from its neigh- 
bors. The farther an agent is away from its nearest neigh- 
bor, the higher speed it has. This strategy can completely 
avoid collisions, however, it results in an aggregated state 
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where the agents move very slow in average. Therefore, 
we further introduce a repulsion that can break down 
the aggregation of agents, and thus sharply speeds up 
the global convergence and enlarges the average distance 
among agents. 



II. MODEL WITH ADAPTIVE SPEED 

We consider each agent as an inelastic ball with ra- 
dius a, limited in a square shaped cell of linear size L 
with periodic boundary conditions. Initially, each agent 
is randomly distributed in the square, with moving direc- 
tion randomly distributed in [— 7r,7r). At each time step, 
the position of the ith agent is updated as: 



X,{t + 1) = + Vt(t), 

and its direction is updated as: 

9,{t + i) = {e,it))r + Ae,, 



(1) 



(2) 



where A6i denotes the thermal noise which is a random 
number uniformly distributed in the interval [— ry, 77] (In 
the main context, we only consider the noise- free case, 
namely ij = 0. A brief discussion about the effect of 
noise is presented in the last section), {9i{t))r denotes 
the average direction of the agents within the horizon 
radius r of the ith agent (including the ith agent itself), 
which reads: 

ta.n[{di{t))r] = {vi sin 9i{t)) r /{vi cos 9i{t))r. (3) 
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In natural swarms, the speed of each agent is alter- 
able, that is, agent may adjust not only its moving di- 
rection, but also its absolute velocity. In the common 
sense, to avoid collisions with other agents, an agent in 
a high-density group should adopt lower speed. Taking 
urban traffic as an example, the speed of an automobile 
is very low in the near-jammed situation, whereas it is 
generally of high speed when sparse automobiles taking 
up the road. Accordingly, we set the speed of the ith 
agent not more than Vi = {di — 2a)/2, where di is the ge- 
ographical distance between two centers of the ith agent 
and its nearest neighbor (see the illustration shown in 
Fig. 1). No matter how the ith agent and its nearest 
neighbor choose their directions in the next time step, 
the restriction can guarantee the distance between them 
no less than 2a and therefore avoid collision. In fact, this 
restriction is not only sufficient, but actually necessary. 
As the direction 9i{t + 1) of each agent in the next time 
step is determined by the average direction within its own 
horizon radius, it is impossible for an agent to know all 
the information of its neighbors, especially the moving 
directions of its neighbors in the next time step. There- 
fore, in order to avoid collision, it is obliged to take into 
account the worst circumstance, that is, two neighbors 
mutually approach. In this worst case, keeping the speed 



(b) 



FIG. 2: Illustrations of locations and velocities in the initial 
configuration (a), and at the 500th time step (b). The param- 
eters are set as L—5, A'^=300, r=l, Vmax=0-03 and a=0.01. 
The length and direction of an arrow represent the absolute 
value and direction of the corresponding agent's velocity. 



of each agent (labeled by i) no more than {di — 2a) jl is 
the only way guaranteeing the absence of collisions. 

Accordingly, the absolute velocity of each agent is up- 
dated with the following rule: 



Vi{t + 1) = Min 



di — 2a 



(4) 



Clearly, when the distance between an agent and its near- 
est neighbor is longer than 2vmax + '^0:, its following speed 
can achieve the maximum; otherwise, its speed is limited 
as {di — 2a) /2. 

Moreover, in order to quantify the consensus of moving 
directions, an order parameter [l8| is introduced as: 



Va 



, < K < 1, 



(5) 
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Especially when Vb = 0, all agents share the same speed. 




FIG. 3: (Color online) Illustration of the motion protocol with 
repulsive effect. 




FIG. 4: The distribution of positions and velocities at 500th 
time step in the scattering model. The parameters are set 
as L=5, Af=300, r=l, Vmax=0.03 and a=0.01. The length 
and direction of an arrow represent the absolute value and 
direction of the corresponding agent's velocity. 

where Vi = \vi\. A larger value of Va indicates better 
consensus. Since the speed in this model is no longer 
constant, it is necessary to introduce another order pa- 
rameter Vb to evaluate the consensus of the absolute ve- 
locity, as: 

Vb = ^^, Vb>0, (6) 

V 

where v = (vi) is the average absolute velocity of all the 
agents, and Aw^ is the variance of the absolute velocity. 
Apparently, a smaller Vb corresponds to better consensus. 



0.7 - 




150 300 450 600 750 900 

b Time Step 



0.032 — ' 




300 600 900 1200 1500 1800 

C Time Step 



FIG. 5: (Color online) Va, Vb, and the average speed versus 
time steps under different repulsion strengthes. The parame- 
ters are set as L=5, A'^=300, r=l, Vmax ~ 0.03 and a=0.0I. 
All the data come form the average results of 500 independent 
runs. 

Numerical results reveal that after the direction con- 
sensus, speed still varies. Figures 2(a) and 2(b) respec- 
tively illustrate the locations and velocities of all the 
agents in the initial configuration and at the 500th time 
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step. After a certain time period from the beginning, 
the positions of agents are not uniformly distributed and 
an aggregation phenomenon appears (see Fig. 2(b)). 
Therein the average speed in a high-density area is much 
slower than that in a low-density area. This aggregated 
state can be understood as follows: agents in a high- 
density region agglomerate together and mutually move 
in a low speed, thus they can seldom disperse apart. 
Meanwhile they take up the way of their subsequent peers 
whose speed is higher, making the high-density area con- 
gregate more agents, and in turn achieving even higher 
density and slower speed (of course, on the other hand, 
the density is limited by the size of agents, a). More- 
over, the motions of agents are similar to the laminar 
flow in hydromechanics: when Va gets close to 1, each 
agent is moving along a line with the same direction and 
will never diverge from its final track. Thus, different 
layers present various flowing speeds. 

In the current model, the nearest distance among 
agents in high-density areas is very close to 2a, making 
the involved agents move in a very low speed (close to 
0); in the meantime the nearest distances among agents 
in low-density areas are usually more than 2a + 2vmax, 
accordingly the involved agents can achieve a high speed 
(close to v„iax)- Consequently, the absolute velocities of 
all agents in the whole system can be in a high diversity. 
Only in a low-density layer can the agents maintain high- 
speed movement in a comparative long term. As a matter 
of fact, the swarm never gc;t speed consensus even with 
identical direction, as shown later in Fig. 5(b). For the 
purpose of making all the agents achieve the consensus 
with higher speed, it is necessary to introduce a certain 
repulsion to avoid agglomeration. In addition, denoting 
Tij — 2a as the safe distance between the ith and the jth 
agents, where rij denotes the geographical distance be- 
tween the ith and the jth agents. In real applications of 
unmanned air vehicles and auto-robots, the longer safe 
distance is favorable. Therefore we hope a properly de- 
signed moving protocol with repulsion could make the 
safe distance longer. 



III. SCATTERING MODEL 

Based on the strategy with adaptive speed mentioned 
above, in this section, we introduce a repulsion to en- 
large the safe distances among agents. We assume: (1) 
the direction of the repulsion should be along the line 
of two agents, and (2) the magnitude of the repulsion 
should decrease with the increase of distance between 
two agents. Moreover, as long as the distance between 
two agents is over 2vrnax + '2a, no matter how they choose 
their directions and speeds, collision will not occur in the 
following steps. Considering this, the repulsion in our 
model should be a short-distance force and work only 
when the distance is shorter than ro (ro = 2vmax 

+ 2a). 



Accordingly, we set the form of repulsion force as: 

. ^ I u X exp(-^37^) X If < ro 
[ nj > ro 

where u is a free parameter. Since the mass of an agent 
plays no role in the present model, we suppose the repul- 
sive effect (caused by the repulsion) can directly affect 
the velocity vector in the next time step (see Fig. 3 for 
an illustration). 
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FIG. 6: Minimal geographical distance between pairs of agent 

in the stable state of the standard Vicsok model, dmin, versus 
the number of agents, A*'. Each data point is the average 
of 1000 independent runs. The restriction to avoid collisions 
with agent size a = 0.01 corresponds to dmin > 0.02. 

After introducing such a repulsion, the moving direc- 
tion of each agent is determined not only by the av- 
erage direction within its horizon radius, but also by 
the repulsive effect. The synthesis of repulsive effect /j 
{fi = X^^i fji^ determined by Eq. (7)) and the average 
velocity Vi (whose direction and magnitude are respec- 
tively determined by Eq. (3) and Eq. (4)) is set as the 
following moving direction of the agent (see Fig. 3). On 
the other hand, the absolute velocity should also follow 
the Eq. (4). Numerical simulations, as shown in Fig. 
4, indicate that this new protocol can effectively scatter 
the aggregated agents (take Fig. 2(b) as an example for 
comparison). Actually, under this protocol, each agent 
can hold a certain distance (much longer than the system 
mentioned in Section II) with its neighbors, and therefore 
achieves its maximal speed, Vmax- 

We also investigate the effects of repulsion strength by 
adjusting the parameter u. Figure 5(a) shows that the 
convergence of moving direction is not sensitive to the 
repulsion strength. However, a larger value of u corre- 
sponds to a shorter time for the system to achieve the 
consensus of speed, as well as a higher average speed in 
the steady state (see Fig. 5(b) and Fig. 5(c)). Consid- 
ering Eq. (4), the larger average speed actually implies 
that the average distance between agents is longer. Note 
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that, when u = 0, Vb can not approach 0, and the average 
speed is very low. 
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FIG. 7: (Color online) Comparison of order parameter Va 
in the Vicsek model and the scattering model under noisy 
environment. The parameters are set as L=5, A''=300, r—1 
and Vmax~0.03. In scattering model, u=0.01 and a =0.01. All 
the data come form the average over 500 independent runs. 
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FIG. 8: (Color online) Comparison on convergence time be- 
tween the standard Vicsek model and the present model (i.e., 
the scattering model) in the absence of noise. In the Vicsek 
model and the present model, the convergence time for Va is 
defined as the required time steps making Va larger than 0.99; 
while the convergence time for VJ, is defined as the required 
time steps making Vb smaller than 10~^. Blue dash curve 
represent the simulation result for the Vicsek model, while 
the black squares and red circles represent the results for Va 
and Vb, respectively. The parameters are set as L—5, r—1 
and Vmax~0.03. In the present model, m=0.1 and a =0.01. 
All the data come form the average over 5 x 10'' independent 
runs. 



IV. CONCLUSION AND DISCUSSION 

As long as we consider the sizes of agents, it is not 
only possible but actually necessary to propose a proto- 
col to avoid collisions among them. Although the Vicsek 
model [isj has achieved a great success in mimicking the 
self-driven swarm, it cannot guarantee the absence of col- 
lisions. We report in Fig. 6 a simple simulation of the 
noise-free Vicsek model neglecting the sizes of agents. 
As the growing of the population, in the stable state, the 
minimal geographical distance between pairs of agents 
decreases quickly. If the size of agent is set as a = 0.01, 
then the minimal distance to avoid collisions must be 
larger than 2a = 0.02. That is to say, the standard Vic- 
sek model can only hold less than 200 agents with size 
0.01 in an 5 X 5 square. In comparison, the current model 
with adaptive speed can hold thousands of such agents. 

However, numerical simulations showed an aggrega- 
tion phenomenon in the current model, which impedes 
the convergence of speed. To overcome this blemish, we 
introduce a repulsion to scatter the aggregated agents. 
The simulation results are exciting: Each agent can hole 
a certain personal space; what is more, they can quickly 
achieve speed consensus and move in a very high speed. 
Numerical results also indicate that the stronger the re- 
pulsive effect is, the less convergent time it takes to 
achieve the consensus. In section II, we have already 
proved that even two neighbors mutually approach, the 
adaptive strategy can still avoid possible collision. There- 
fore, in any event, collision will never occur in the scat- 
tering model. 

Furthermore, it is well known that the thermal noise 
can also play a significant role in determining the moving 
directions of agents. Thus, we need to check whether our 
rule is robust in the presence of noise. The numerical 
result indicates that even in the noisy environment, in 
the stable state, the average distance and average speed 
are both larger than those without the repulsion. The or- 
der parameter for direction consensus of course decreases 
with the increasing of noise strength, 77, and it exhibits 
almost the same trend as the standard Vicsek model (ac- 
tually, it is a little bit larger than the Vicsek model, see 
please the simulation result shown in Fig. 7). 

In the noise-free Vicsek model, given r and L, the 
convergence is faster with more agents (i.e., larger N) 
since they will have more frequent communications in a 
denser circumstance. Actually, a recent numerical study 
[2^ indicates that the convergence time scales approxi- 
mately as (InA^)"^ '^, that is, the larger the population 
is, the shorter the convergence time is. In Fig. 8, we re- 
port the simulation result on the convergence time in the 
noise- free Vicsek model (see the blue dash curve) , where 
the threshold quantile is set as Va = 0.99. It decreases 
monotonously with the increasing of N, in accordance 
with Ref. [28|. In contrast, in the present scattering 
model, more effort should be taken to avoid collisions 
in the denser circumstance. Figure 8 compares the con- 
vergence time between the standard Vicsek model and 
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the scattering model in the absence of noise. One can 
find that in the sparse circumstance, N < 600, the con- 
vergence times of the Vicsek model and the scattering 
model are almost the same, while in the denser range, 
the convergence time in the scattering model quickly in- 
creases versus the slowly decreasing of that in the Vic- 
sek model. The convergence time for absolute velocity 
increases even most quickly than that for moving direc- 
tion. This result indicates a limitation of the present 
model, namely it can not efficiently deal with the sys- 
tems with huge population. Accordingly, how to design 
an efficient method to simultaneously guarantee the ab- 
sence of collisions and the quick convergence is still an 
open problem for us. Anyway, in the case of a = 0.01, 
the standard Vicsek model can avoid the collisions only 
if the number of agents is less than or about 100 (see 
Fig. 6), while the scattering model can hold about 600 
agents with the same speed of convergence. We therefore 
believe the scattering model can find applications in the 
design of motion protocol for self-driven agents. 

Some difficult yet important problems about the con- 
servative model remain to be further explored. For ex- 
ample, if the ahead ones of a group of agents need not 
pay attention to the following ones (that is, each agent 
only receives information in a sector ahead in its moving 
direction rather than all the neighbors within its sight 
radius [llj), collisions may automatically disappear. If 
the swarm needs shorter time to get convergence while 
avoiding the collisions, it may indicate that the complete 



communication is not always the most efficient manner 
while partial communication may be better in some cases 
[29I . SH . In addition, the properties of the phase transi- 
tion induced by the noise (see, for example, in Ref. [33|, 
Gregoire and Chate showed that a swarm model with 
repulsion as well as the minimal Vicsek model suffers 
a discontinuous phase transition) remains an open issue. 
Though not the focus in this article, it worths an detailed 
investigation in the future. 
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